Conditions for extreme sensitivity of protein diffusion in 
membranes to cell environments 



Yaroslav Tserkovnyak 
Lyman Laboratory of Physics, Harvard University, 
Cambridge, Massachusetts 02138, USA and 
Department of Physics and Astronomy, 
University of California, Los Angeles, California 90095, USA 

David R. Nelson 
Lyman Laboratory of Physics, Harvard University, 
Cambridge, Massachusetts 02138, USA 
(Dated: February 6, 2008) 

Abstract 

We study protein diffusion in multicomponent lipid membranes close to a rigid substrate sepa- 
rated by a layer of viscous fluid. The large-distance, long-time asymptotics for Brownian motion 
are calculated using a nonlinear stochastic Navier-Stokes equation including the effect of friction 
with the substrate. The advective nonlinearity, neglected in previous treatments, gives only a small 
correction to the renormalized viscosity and diffusion coefficient at room temperature. We find, 
however, that in realistic multicomponent lipid mixtures, close to a critical point for phase sepa- 
ration, protein diffusion acquires a strong power-law dependence on temperature and the distance 
to the substrate H, making it much more sensitive to cell environment, unlike the logarithmic 
dependence on H and very small thermal correction away from the critical point. 
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We revisit here the problem of particle diffusion and relaxation of concentration fluctua- 
tions within a two-dimensional incompressible hydrodynamic medium in proximity to 
a rigid substrate . Interest in the two-dimensional Brownian motion has been largely 
motivated by protein diffusion in cell lipid membranes, with essential features captured by 
the seminal theory of Saffman and Delbriick 
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who predicted logarithmic dependence of 
;he diffusion coefficient on a large-distance cutoff. The central result of the linearized theory 
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j'or the translational diffusion coefficient is 

D~^'-, (1) 

assuming an appropriate large-distance cutoff I* 3> a, where a is the small lateral size 
of the diffusing protein. See Fig. [T] for a sketch of the physical situation. Here, ks is 
the Boltzmann's constant, T the absolute temperature, h the membrane thickness and [i its 
dynamic viscosity. The cutoff I* depends on membrane dimensions and also on environmental 
factors such as the viscosity of surrounding fluid or the distance to a nearby rigid wall. Eq. 
allows an estimate of the protein size in terms of a known viscosity or vice versa, if one can 
selectively probe diffusion of a givenprotein. Protein diffusion coefficients were measured, 
for example, by Peters and Cherry p in large lipid vesicles by monitoring fluorescence of 
diffusing proteins after a circular area of the membrane was illuminated by a laser. Owing 
to the logarithmic dependence of the diffusion coefficient on a, however, the precise form 
of the diffusion coefficient in terms of basic parameters of the problem is clearly required in 
order to make such techniques useful. Possible I* (and thus D) sensitivity to proximity to 



a rigid substrate can be used to investigate t 
to a nearby substrate such solid wall 



re physics and biology of membrane coupling 
Proximity to a substrate affects transport 



and possibly also structural properties of the membrane, effects that could be biologically 
relevant for cell signaling and sensing of the external environment Q. 

For an isolated flat two-dimensional membrane with lateral momentum conservation, I* 
is set simply by the membrane size. In practice, lateral momentum conservation requires a 
negligible viscosity of the surrounding fluid and an appreciable offset from any rigid walls. 
When I* — > 00, the diffusion coefficient Q diverges, which is reflected in the logarithmic 
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factor in the mean-square displacement of a diffusing particle: (r 2 (t)) cx tint |ll,l2fl- I n this 
case, the diffusion coefficient becomes effectively time dependent, growing logarithmically in 
time and eventually saturating at the value given by Eq. (JTJ) with I* determined by geometric 
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substrate 

FIG. 1: Schematic of a flat two-dimensional lipid membrane of thickness h separated by a layer 
of a bulk liquid of depth H from a parallel rigid substrate. Since realistic lipid bilayers do not 
permit shear across the membrane, the velocity field in the membrane is taken to be strictly two- 
dimensional, characterized by a (two-dimensional) dynamic viscosity fj,h. The liquid is described 
by a dynamic viscosity ji' (for water, // <C fj). No-slip boundary conditions are assumed at both 
the substrate and membrane interface. Two scenarios for thermally-driven protein diffusion are 
considered in the text: We first discuss Brownian motion of proteins of lateral size a within a uni- 
form lipid background. In the second model, proteins segregate in lipid "rafts" of (predominantly) 
a particular chemical composition above the critical temperature T c of a binary lipid mixture. The 
large-scale diffusion behavior of proteins is then reformulated in terms of the raft diffusion of size 
£. Physically, £ is the correlation length of the binary mixture, which diverges as a power law 
near the critical temperature. As £ grows near T c , protein diffusion becomes extremely sensitive 
to temperature and substrate depth H, see Eq. (|23j). 

constraints (either a finite size or radius of curvature). The divergence of the linear mobility 
when I* — > oo is the Stokes paradox There is no solution of the low Reynold's number 
linearized-flow equations for steady translational motion in two dimensions. The time- 
dependent diffusion coefficient D and mobility b are however still related by the Einstein 
relation: 

D = bk B T . (2) 
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Upon taking into account the inertia of the viscous fluid (in the form of the advection), the 



a. 



nonlinear mobility becomes finite but logarithmically growing with the inverse velocity 
A cutoff I* for the large-distance divergence of the diffusion coefficient and corresponding 
linear mobility related by Eq. (j2J) can be also provided by attaching the membrane to a 
momentum sink on either or both sides. Embedding the membrane into a slightly viscous 
three-dimensional liquid with viscosity // <C //, for example, provides a cutoff /* — > R* = 

HQ. 

For biological membranes with \ij ' jj ~ 100, this is shorter than geometric 
cutoffs, provided we take h ~ 10 nm and consider cells larger than /im. 

Evans and Sackmann jj| discussed the problem of a coupling to a nearby rigid substrate, 
see Fig. [TJ (Such a setup was used in Ref. P to build a two-dimensional shear viscometer.) 
When the coupling is mediated by a thin layer of viscous fluid of depth H <C R* and viscosity 
//, the diffusion coefficient becomes 



D k ° T 



a K\(a/5) a 2 n 1 
5K (a/5) + 25* 



(3) 



where K and K\ are modified Bessel functions of the second kind of orders zero and one, 



respectively, and 5 = V R*H. We plot Eq. (JHJ) in Fig. |2j When 5 ^> a, Eq. (JHJ reduces 



evant cutoff 5(H); 
and the cutoff /* 
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5, 



to the form (0) with /* = 5. Increasing H naturally increases the re 
however, the assumptions leading to Eq. (JHJ) break down when H > R 
saturates at 5(H = R*) = R*, where the remote rigid substrate becomes unimportant. In 
the opposite limit of 5 < a (while still H < R*), Eq. © reduces to D w k B TH/(27ra 2 fi'), 
from which one easily recognizes the mobility determined by the direct laminar-flow friction 
of theparticle with the substrate. In these and similar linearized low-temperature theories 
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flow equations are solved for specific geometries. For example, the diffusing particle 
can be assumed to be a rigid cylinder of height h and radius a (with, e.g., a no-slip or 
no-tangential-stress boundary condition for liquid flow In the following, we present 

an approach that focuses on the long-time large- distance physics from the onset, without 
reference to any short-distance details, apart from a characteristic protein size a. 

Our interest in the problem is twofold. On the one hand, we want to understand the inter- 
play between thermal fluctuations and advective nonlinearity omitted in previous treatments 
of protein diffusion and linear mobility. In addition, we wish to investigate the behavior 
near a critical temperature T c of the multicomponent lipid mixtures that make up biological 
membranes. One can efficiently deal with nonlinearities using the dynamic renormalization- 




FIG. 2: The solid line shows the undimensionalized diffusion coefficient (J3J) in units of ksT / '{Air^h) 
as a function of 8/a. The dashed (dotted) line represents the large (small) 8/a limit discussed in 
the text. 



group (RG) approach jlO, 11, 12, 13], which will allow us to scale the large- distance long-time 
physics onto the "asymptotically free" (i.e., linear) theory. We will then apply some general 
ideas from the theory of dynamic critical phenomena [12| to the lipid and protein diffusion 
close to a critical point. Of a central interest to this paper is the sensitivity of these effects 
to proximity of a rigid wall, which exerts friction and bleeds momentum from the mem- 
brane. We thus want to generalize the result ® to include the interplay between nonlinear 
corrections and thermal noise and to also adapt the theory to include critical fluctuations 
in multicomponent lipid membranes. 

Before presenting details of the model and calculation, let us first outline what we can 
anticipate based on some simple but general arguments 1^. Suppose the binary lipid 
membrane is close to a critical point for phase separation. Below the critical temperature T c , 
lipids form laterally-phase-separated regions ( "rafts" ) . Specialized proteins can segregate in 
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lipid rafts for cell signaling, protein organization, and endocytosis. Raft-like structures were 
also observed in artificial lipid membranes consisting of a binary mixture with one phase rich 
in the unsaturated lipids and the other rich in the saturated lipids and cholesterol jl4. 15, 



Below T c , phase-separated membranes also have inhomogeneous elastic properties, which 
lead to formation of nontrivial spatial structures [h! Q| that are beyond our analysis. 
For simplicity, we focus on membrane diffusion at temperatures slightly above T c , where 
lipids form nearly uniform mixture with correlation length diverging on approaching T c . If 
proteins strongly prefer one kind of lipid versus the other, they will preferentially segregate 
in the correlated regions (which we can also consider as "rafts") corresponding to one of 
the emerging phases. On time scales longer than it takes for the protein to diffuse across 
such rafts, the overall protein motion becomes essentially passive and governed by lipid 
concentration intermixing. It will be shown that the protein size a must then be replaced 
by the correlation length £ in order to obtain the effective diffusion coefficient. The latter 
is then approximately given by the Einstein relation (j2j) with mobility b scaling as £ 2_d //i 
(assuming momentum conservation within the membrane), which is inversely proportional 
to the drag force exerted by the membrane fluid on a drifting raft of size £, where we in 
general consider a ^-dimensional system (eventually setting d — 2). We thus have 

Dfi oc k B TZ 2 - d , (4) 

where in general both D and \x can be temperature dependent. For d > 2, the right-hand 
side of Eq. (£Q) vanishes near the critical temperature where £ diverges. Short-range physics 
on length scales shorter than £, which is disregarded here, then has to be invoked in order to 
calculate the product Dfi. When d < 2, on the other hand, Eq. (J3J diverges near T c signaling 
emergence of important long-distance processes, d — 2 is the critical dimension where one 
may expect only logarithmic singularities due to long-distance physics, as we will indeed see 
in the following. When there is a rigid substrate exerting a friction on the two-dimensional 
membrane dynamics, however, the diffusion coefficient can acquire a power-law dependence 
on both £ and the distance to the substrate, making the diffusion very sensitive to the 
environment near T c . In particular, as explained in the following, there is a regime where 
the protein size a is effectively replaced by the correlation length £ in Eq. @, resulting in 
Eq. (jZH|) for large £. Note that in this case D depends on the viscosity of the surrounding 
fluid and not the membrane viscosity itself. 
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Eq. Q is also relevant for diffusion away from the critical point, where £ may be replaced 
by 1/q and the diffusion coefficient and viscosity become wavevector-dependent functions, 
D(q) and p(q). The long-distance physics at q — > is important for d < 2, consistent with 
the result of the linearized theory for d = 2, which gives 

Dp oc ln(l/ga) (5) 

for wavelengths shorter than I* [obtained by replacing I* with 1/q in Eq. JJ}]. We will show 
that although nonlinearities and thermal fluctuations leave Eq. (JHJ) unmodified away from 
the critical point, the logarithmic singularity is not entirely contained in D; in fact, the 
effective viscosity p may also acquire a finite-temperature logarithmic q dependence. We 
study this phenomenon first, before discussing diffusion near the critical point, in order to 
understand how the membrane viscosity itself depends on temperature and relevant length 
scales. 

Let us now formulate the model. Focusing on the long-wavelength two-dimensional dy- 
namics, we write Navier-Stokes equation for membrane fluid velocity field subject to a ran- 
dom forcing function f(r, t): 

d t v + A(v • V)v = ~Vp - av + z/V 2 v + f , (6) 
P 

where all vectors have two components in the xy plane, e.g., r = (x,y). The field v(r, t) is 
the position- and time-dependent fluid velocity, p is a constant mass density and v = p/p the 
kinematic viscosity of the membrane fluid. (In writing a strictly two-dimensional Navier- 
Stokes equation, we have assumed a large planar shear viscosity within the membrane, 
relative to the fluid outside.) The coupling a describes friction due to a nearby rigid sub- 
strate. If the friction is mediated by a constant-gradient laminar flow of a liquid layer of 
thickness H with viscosity //, then a = p'/(hHp). (For simplicity, we introduce the cou- 
pling a only on one side of the membrane.) We are assuming here that R* = (p/p')h^> H; 
otherwise the membrane does not sense the substrate and H has to be replaced by R* in 
the definition of a for lateral lengthscales smaller than R*. The pressure term — (l/p)Vp 
enforces incompressible fluid flow: 

V ■ v = . (7) 



The second term on the l.h.s. of Eq. © is the usua 



We will be able to treat it using an RG approach jll] where A becomes a perturbative param 



(nonlinear) advective term when A = 1. 



eter weakly irrelevant for large-distance asymptotics. The fluctuation-dissipation theorem 
dictates the form of f in correspondence with the dissipative terms in Eq. ©: 

(/i(k, u})f s {\t, u/)) = 2r(a + z/A; 2 )P lJ (k)(27r) 3 5(k + k')S(ou + u') , (8) 

where f (k, to) = J d 2 rdti(r, t) exp(iut — ik • r) is the Fourier transform of the velocity distri- 
bution, r = k B T I {hp) is a rescaled temperature, and -Pij(k) = 5^ — kikj/k 2 is the transverse 
projection operator. We note that when coupling a is mediated by a thin liquid interlayer 
connecting the membrane with the substrate, we account for it by assuming a simple laminar 
flow, which is not perturbed by thermal fluctuations and which responds instantaneously to 
the driving shear provided by the membrane dynamics. The nonlinearities and stochasticity 
are fully contained within the membrane. We thus neglect quasi-two-dimensional nonlin- 
ear and thermal effects as well as inertia in the interlayer. However, we expect that these 
complications simply lead to a renormalization of parameters such as A, p, a, and v in 
Eq. ©. 

The particle (protein) concentration c is treated in the dilute limit as a passive scalar 
characterized by some local diffusion coefficient Do and statistically advecting along fluid 
flows (with no back- reaction on the hydrodynamics) described by Eqs. (ffi|)-([5]): 

d t c + A(v- V)c = A)V 2 c, (9) 



which can be thought of as the Fokker-Planck equation for the probability distribution of 

At this point, we disregard protein 



an individual particle undergoing Brownian motion ^ | . 
coupling to the rigid substrate and assume that it is felt only indirectly through the fluid 
friction a and particle's drag by the fluid flow. Therefore, A = 1 in Eq. Q. As explained 
in the Supplemental Material, this coupling then obeys the same RG flows upon coarse 
graining as A in Eq. ©: that is why we are using the same variable A in the two cases. 
For sufficiently shallow depths H, however, the fluid-mediated particle-substrate friction 
becomes important, as already mentioned in the linearized regime described by Eq. (J3j). 
This can be easily taken into account in the final result for the mobility [and thus diffusion 
coefficient, according to the Einstein relation (J2J)], as an additional drag by the rigid substrate 
on the diffusing particle. In the following, we solve Eqs. focusing on the viscosity 

and diffusion-constant renormalization in the infrared limit. 

Let us first calculate viscosity renormalization due to the advection term in Eq. ©. 
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Proceeding along the lines of Ref. Ill ], we Fourier transform Eq. (JUJ) to obtain 

i f cPqdfl 

vi(k,u) = G (k,u)fi(k,u) - -XP Hj (k)Go(k,uj) J ^(q, Sl)vj(k - q, w - ft) , (10) 

where Puj(k) = kiPji(k) + kjPaik), G (k,uj) = (—iuj + a + vk 2 )' 1 is the response function 
of the linear theory with A = 0, and summation over repeated indices is implied. Eq. ()10|) 
can be solved perturbatively in A by iteration. We first introduce a short-distance cutoff 
1/A for the velocity field v and the stochastic field f. Physically, 1/A is of the order of 
the interlipid separation, the scale at which the continuous hydrodynamics obviously fails. 
Correspondingly, the Fourier transforms must be cut off at momentum A. The lowest-order 
correction to the full response function G(k, u) due to the advection term oc A results in the 
following contribution to the wave- vector-dependent viscosity (setting A = 1) [ll| : 

Av » — 1— ln(A/fc) (11) 
Wilis 

at uj = 0, assuming k is larger than the momentum cutoff 1/5 = yvja. = a/ (Jl] Jj/JhH set 
by a finite a. When k <C 1/5, the argument of the logarithm in Eq. (fTTj) must be replaced 
by A5. The importance of advection due to thermal fluctuations described by Eq. (|TT| thus 
hinges on the size of the prefactor r/(167rz/). Taking room temperature, water density, and 
h ~ 10 nm for the evaluation of r = ksT/(hp) and setting membrane viscosity v to be 
100 times larger than that of water, we get t/(8ttu 2 ) ~ 10 -9 . Clearly this gives only a tiny 
correction to the linearized theory for biologically-relevant problems. 

It is nevertheless interesting to pursue a nonperturbative calculation in the infrared limit 
k — > for a very small damping a — > leading to a logarithmic divergence of the correction 
(fTTj) . This exercise allows additional insights into the peculiarities of the two-dimensional 
dynamics. In order to discuss the nonperturbative asymptotics, we proceed using the dy- 



namic RG 



131 ] in two standard steps: First, we reduce the cutoff to Ae~ l , where / 1, 
and rescale parameters a, v, and A of the Navier-Stokes equation (@). The requirement 
is that going from initial (physical) set of parameters (A,a,^, A) to [Ae~ l ,a(l),u(l), X(l)] 
does not modify velocity correlators and response functions averaged over the random force 
(jEJ)- Second, we rescale momentum by e l , so that the cutoff is unchanged by the two steps, 
and also conveniently choose frequency rescaling. The velocity field is rescaled by requiring 
that the relation (jHJ) remains unmodified in terms of the /-dependent a and v. Finally, the 
parameters a, z/, and A are correspondingly renormalized. This procedure is then repeated 



continuously with accumulated rescaling exponent I going from to a chosen value much 
greater than 1. 

Since we adapted the scheme developed in Refs. 11, 3> we wm av °id repeating the details 



of the derivation. A summary is given in the Supplemental Material. Solving Eq. (J31|) in 
the Supplemental Material gives for the renormalized viscosity at k — »■ 



z/2 + ^ln(A5 r ), (12) 

07T 



where 



S r =^hH (13) 

itself depends on the renormalized (dynamic) viscosity fi r = v T p. If the membrane is nearly 
isolated, i.e., a = f//(hHp) « 0, the second summand under the square root in Eq. (JT2Jl 
dominates and the renormalized viscosity grows logarithmically when k — > 0, eventually 
saturating at the value determined by the distance cutoff S r : 



Vr(k^0)*J-\nl-)^J-\n(A6 r ), (14) 




which agrees with the result obtained in Ref. [llj . In practice, however, as we have dis- 
cussed, this limit is hard to achieve. The thermal renormalization of the viscosity due to 
the advective term in the Navier-Stokes equation (JHJ) turns out to be very small and instead 
v r ~ v. 

Having studied the large-scale long-time velocity-field dynamics, we will now discuss con- 
sequences for the passive-scalar advection described by Eq. (JHJ). This equation can be solved 
iteratively treating A as a perturbative parameter, similarly to the velocity-field expansion 
(I10J) but with a slightly different physical formulation: For the velocity field v(r, t), we 
were interested in the response to thermal and external forces driving the fluid dynam- 
ics; for the particle diffusion, on the other hand, we do not have any particle sources or 
sinks, but instead want to consider the evolution of the initial distribution c(r, t = 0) at 
t > 0. It is thus natural to Fourier transform c(r,t) in space but Laplace transform in 
time: c(k, uj) = J d 2 r J °° dt exp(iut — ik- r), where u has an infinitesimal positive imaginary 
component. The scalar field c(k, uj) is then given according to Eq. Q by 

c(k, u) = G (k, uj)c(k, t = 0) - i\G (k, u)ki J T^F c ( q ' _ °L> u ~ fi ) > ( 15 ) 
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where Qo(k,u) = (—iu + Dok 2 )^ 1 is the diffusion propagator [v(k, uj) is still the full Fourier 
transform, as before] . The technical details for solving Eq. (|T3j) are once again summarized in 
the Supplemental Material. Interestingly, we find that D T ~ ln 1 ^ 2 (5 r /a) at the very longest 
lengths due to thermal fluctuations, in contrast to Eq. However, at the temperatures 
of interest for membranes, we find in practice that v r pa v and v r ^> D r . Eq. (|34j) of the 
Supplemental Material then gives for k — > simply 

U r «JLlB^V (16) 



Attu \ a t 

neglecting the short- distance contribution D and assuming 5 r 3> a. This is exactly the 
Saffman-Delbruck result (Q) for the relevant distance cutoff /* = 5 r . We furthermore make 
the ansatz that at long times and distances for the diffusion problem one can in general use 
the results of the linearized theory, such as Eqs. (0) and Q, as long as the renormalized 
viscosity v r and self-consistently-found distance cutoff 5 r [Eq. (JT3J)] are used. We expect this 
to be correct, up to numeric factors of order unity (such as the factor 2/ A in the Supplemental 
Material), even in the extreme high-temperature limit when thermal fluctuations dominate 
equilibrium properties. The final result for the renormalized diffusion constant in the system 
sketched in Fig. ^is thus a trivial generalization of Eq. (jH)): 



a Ki(a/5 r ) a 



2 



-1 



(17) 



Anfi r h l5 r K (a/5 r ) 25% 

Like Eq. ©, this result includes the fluid-mediated friction between protein and the sub- 
strate, which we have not explicitly considered in the above discussion. As we have already 
explained, we can neglect thermal renormalization of v T (and thus 5 r and fi r ) for biological 
systems. 

We now discuss how these results are modified close to a critical point for phase separa- 
tion in binary lipid membranes, leading to a much more dramatic dependence of diffusion 
on temperature and membrane environment. To be specific, let us consider concentration 
interdiffusion in a s ym metric binary mixture of, say, lipids A and B with concentrations 
c^(r, t) and ce(r, t) \12\ . By symmetry, the critical point corresponds to equal average con- 
centrations: (ca) = (cb) = c/2, where c is the averaged total concentration (taken to be 
constant). The order parameter ip(r,t) is given by the normalized difference of the concen- 
trations: ipijc^t) = [ca{t, t) — cs(r, t)]/c. The total concentration, c^(r, t) + ce(r, t), does 
not couple to the ip diffusion mode in this case. Following Ref . |l2| , we introduce a fictitious 
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external potential which couples to ip, by choosing opposite potentials for the two lipid 
components: \ia = $ and fi B = — The force on a region with a roughly uniform order 
parameter ip of size £ (the order-parameter correlation length) is then given by 

f 4 = - VF = -£ 2 ^V$ , (18) 

assuming a long- wavelength (on the scale of £) driving field $, where F is the free energy. 
The region will accelerate until reaching velocity v when f<j, is canceled by the viscous drag 
f„ = — v/b characterized by mobility h: £ 2, 0V<& + v/b = 0. We thus get for the net current 
density of the order parameter ip within the region of size £: 

= ^ v = -6£ 2 ^ 2 V$ = -aV$ , (19) 

where 

a = b^tfj 2 = bk B T X (20) 

is the concentration conductivity associated with the potential $ and x is the susceptibility 
of the concentration difference if). The last equality in Eq. (|20j) is obtained by invoking the 
fluctuation-dissipation theorem, 

(F) = (k B T/e)x, (21) 

implying averaging over correlated regions (over the entire system). For a smooth static 
external potential, ip = in equilibrium, and thus the current density must be pro- 
portional to V($ + t/>/x). In the absence of the fictitious potential $, the concentration 
flows are then described by the diffusion coefficient D: = —D'Vip with 

D = — = bk B T , (22) 

X 

where we have used Eq. (j2~U|) for the last equality. Since ip is a conserved quantity, d t ip = 
— ^ 'Ji> — DV 2 ip. Eq. (|2*2*|) is simply the Einstein relation (j2J), with one conceptual subtlety: 
D and 6 describe not the transport coefficients of individual small particles, but rather 
the (inter) diffusion coefficient of a binary mixture and the mobility of a region of size £, 
respectively. In particular, we can approximate the latter by the mobility of a "particle" of 
size £, which in turn corresponds to the diffusion coefficient p7jl . We conclude from this line 
of argument that the diffusion coefficient of proteins "riding" rafts with order parameter ip 
of a given sign [determined by protein affinity to a given kind of lipids, either A (-0 > 0) 

12 



or B (?p < 0)] is given by equation (JTTj) with a — > £. In particular, as T — > T c and the 
correlation length £ diverges, £ 3> S r and the diffusion coefficient approaches 



k B TH 
2tt£V 



(assuming if <C i?*), as we discussed in the introduction (see also the small-<5 asymptotic 
of Fig. EJ). The membrane viscosity fi in fact also has a diverging contribution near T c , but 
fortunately with a very small exponent (~ 0.16 in two dimensions) [12], so that as long as 
£ 5 r it drops out of Eq. (|2~3J) for the diffusion coefficient. The exponents that govern the 
correlation-length and susceptibility 

X (T) oc £(T) 2 ^ (24) 

divergence near the critical temperature can be calculated in the two-dimensional Ising 
model (which is in the same universality class as the binary mixture), giving £ oc (T — T c ) _1 
and rj = 1/4. Upon combining Eqs. (j22j) . (J2*3*j) . and (J21J), we find for the conductivity scaling 
cr(T) oc C,(T)^ V , which vanishes with a small exponent 77 at the critical temperature. In 
practice, however, equation (|2T)|) with b standing for the mobility of a rigid region of size £ 
neglects short-range (on the scale of £) fluctuations that can lead to a finite conductivity 
o"o a t the critical temperature. The diffusion coefficient (|2*2*|) then has to be correspondingly 
corrected, giving near T c (with £ ^> 5 r ) 



k B TH a fa \ 2 v 



with the second contribution eventually dominating very near T c . The constants <To and %o 
are nonuniversal, dependent on nonuniversal details of the lipid mixtures; ao is a character- 
istic lipid spacing. 

We conclude with a discussion of one final point: Although the correlation length grows, 
density fluctuations eventually disappear at the critical temperature and the protein diffusion 
is not critically slowed down. Upon combining Eqs. (j2~Tj) and (J2H), we S e t (4> 2 ) — Cf^o/O^ 
near T c , in terms of some dimensionless constant C. We make a mean-field approximation 
for the free-energy difference for a single protein residing in region A (ip > 0) and region B 
(ip < 0), SF ~ K,5ip/2, where 5ip is the order-parameter difference between the regions, and 
assume that the proteins will thermally hop between the correlated regions, and thus freely 
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diffuse within the membrane when k B T ~ k5i/j/2 ~ KyC(ao/^) r1 ' 2 , i.e., when the correlation 
length becomes larger than 

£ max ~ a , (26) 

a relation which in two dimensions is characterized by a large exponent 2/rj = 8. Very close 
to the critical point where the correlation length is larger than Eq. 1)26)1. the protein diffusion 
is thus not profoundly affected by the critical behavior apart from the weak singularity in 
the membrane viscosity. 

We can make a crude estimate of £ max as follows. Suppose k is of order the energy 
difference between a typical protein embedded in the two-chain lipids and a typical protein 
embedded in the one-chain lipids. A typical protein size is 50 A; with a lipid head size of 
6-10 A, maybe 20 lipids (ten from each leaflet) surround each protein. Assume the protein 
spans the bilayer, so each of ~ 12 CH 2 groups of the lipid chains above and below can interact 
with the protein. A typical hydrogen bond or van der Waals attractive interaction with a 
CH 2 might be 0.2k B T at room temperature. Since roughly 100 CH 2 groups are in contact 
in one case (single hydrocarbon chain), and 200 CH 2 groups in the other (two hydrocarbon 
chains), the energy difference is k ~ 20k B T. In situations where both lipid species have two 
hydrocarbon chains, the energy difference could arise from the different head groups (not 
included in the above) and from the typical fact that one species has an unsaturated bond 
in the hydrocarbon chain (making protein interactions with part of the chain difficult) and 
the other does not. Also, the chain lengths might differ by 2-3 CH 2 groups. All these factors 
together might lead to 

« ~ 10k B T (27) 

in less favorable situations. The constant C can be estimated by the mean-field approxima- 
tion of the Ising model |17]: Using the form L = at^ 2 + bip 4 /2 + -f(Vip) 2 /2 [t = (T — T c ) /T c 
is the reduced temperature] of the Landau free-energy density, one gets for the mean-field 



response function x — Il an d the correlation length £ = ^ r y/2at. According to Eq. 
then C = ksT/'-f and rj = (instead of the correct rj = 1/4 since the mean-field approxi- 
mation captures the exponents correctly only at more than four dimensions). For the Ising 
model, a mean-field estimate gives 

C = k B Tl~i ~ (a /R) 2 ~ 1 (28) 
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for short-range interactions, where R is the range of exchange coupling. Using the estimates 
(J27)l and (pS|). we finally conclude that £ max ~ lO^cto according to Ecj. (]26j) . which clearly 
suggests that the reasoning leading to Eq. (|25|) holds for a large range of the correlation length 
£ (apart from temperatures very close to T c where the critical slowing down is suppressed). 

In summary, we have considered the effects of thermal noise, hydrodynamic nonlinearities, 
and proximity to a critical point on protein diffusion in multicomponent lipid membranes 
close to a rigid substrate. We have found that nonlinear advection in the Navier-Stokes 
equation in practice gives only a small correction to the Evans- Sackmann result (jSJ) for a 
uniform lipid composition. The protein diffusion is then only logarithmically sensitive to 
the distance to the wall, assuming the protein size a <C 5 = a/ (fi/fi')hH. Upon approaching 
the critical temperature for the phase separation in multicomponent membranes, however, 
the long-time asymptotic of protein diffusion can be reformulated in terms of the lipid 
"raft" diffusion of size £, which diverges at T c . The diffusion coefficient then becomes 
linearly sensitive to the distance to the substrate H and has a roughly l/£ 2 scaling with 
the correlation length, near T c (apart from a narrow temperature window very close to T c 
where proteins thermally hop between the rafts). If it were functionally advantageous for 
cell sensing and signaling, one might imagine that lipid composition could evolve to be near 
a critical point. 
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SUPPLEMENTAL MATERIAL 



Diagrams defining the intermediate recursion relations (before momentum rescaling) for 

the propagator and vertex following from Eq. (jl(Jj) are schematically shown in Fig. El The 

resulting renormalizat ion-group (RG) flows after momentum rescaling are given by 

da (I) ,,\dj(l) 
- a{l)- 



dl dl 



du(l) 

dl 
d\(l) 

~dT 



!/(/) 



did) 

dl 

' d-yjl) 
dl 



r A(/) 2 A 2 
167ra(/) + z/(/)A 2 ' 

(29) 



+ 



where 7(Z) is an arbitrary exponent for frequency rescaling. It turns out that only viscosity is 
modified by reducing the momentum cutoff A (in particular, the three vertex-renormalization 
diagrams in Fig. El cancel, as in Ref. |nj which had a = 0). A convenient choice for 7(Z) 
ensures that a (I) + u(l)A 2 = const. It is natural to define the renormalized viscosity by the 
velocity correlator 

G«(k, w) = , K I ' V ' f „ = 2rP«(k)Re ^ . 30 

By the fluctuation-dissipation theorem, this is the same viscosity that enters the full velocity 

propagator. Upon expressing the velocity v, wave vector k, and frequency to in terms of 

the rescaled quantities, we can calculate if the renormalized parameters flow into a 

solvable regime. If u — (which will be assumed in the following), RG flows have to be 

stopped before e l k = A when the rescaled momentum hits the cutoff. When k — > 0, I can 

thus be sent to infinity, and the theory becomes asymptotically free as solving Eqs. (|29|) 

gives X(l — > oo) — ► 0. Within perturbation theory, small remaining A only results in a k 2 

contribution to the viscosity v r , which can be safely disregarded in the large-distance limit. 

We thus have to solve Eqs. (J29|) as I — > ln(A/fc) to find a(l) and which define the linear 

velocity correlator that can be mapped onto the full correlator (J3U|) of the original problem. 

The solution of this problem for the renormalized viscosity can be recast as a solution of the 

following first-order differential equation: 

dv T {k) = t k 

dk 16tt a + u r (k)k 2 ' 1 ' 

with initial condition v r (k — A) = v. The renormalized viscosity u r (k) is given by u r (k = k). 
Solving Eq. (}3*T]) results in Eq. (JT2|) in the limit of k — > 0. 
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FIG. 3: Graphical representation of recursion relations for the propagator G (which determines 
substrate friction a and viscosity v) and vertex proportional to A. Single lines are propagators G 
with momentum cutoff A, open circles are averaged force correlators ((HI) on the momentum shell 
Ae - ' < q < A, and bare vertices are given by — (i/2)\Pnj(k) according to the last term in Eq. (|10|). 
Double lines are propagators with scaled cutoff, to order A 2 , and dashed circle is the renormalized 
vertex, to order A 3 . The intermediate momenta are integrated over the shell Ae~ l < q < A. 

The diagrams representing recursion formulas for the diffusion coefficient D and vertex 
A corresponding to Eq. (fTK|) are sketched in Fig. 0] [ijj]. As before, nonlinear corrections to 
the vertex vanish, and we can appropriately rescale c(k, u) in the second RG step so that A 
flows exactly as in Eqs. (|29|) [that is why we used the same A in Eqs. (JBJ) and Q]. We thus 
have to supplement Eqs. (|29p with only one additional flow equation: 



dD{l) 
dl 



D(l) 



dl{l) 
dl 



+ 



T 



A(/) 2 A 2 



4vra(/) + [u(l)+D(l)]k 2 ' 



(32) 



Upon defining the renormalized diffusion coefficient by 

c(k,i = 0) 



(c(k,w)) f 



f -iuj + D r (k,uj)k 2 
and setting uj = 0, we finally obtain the differential equation similar to Eq. 



(33) 



dD r (k) 



T 



k 



(34) 



dk 4tt a + [u r (k) + D r (k)]k 2 ' 

with initial condition D r (k ~ 1/a) = D (physically meaning that the diffusion- const ant 
flow is "turned on" at momenta smaller than inverse protein size, assuming k 1/a). 
D r (k = k) then gives the renormalized diffusion coefficient D r (k), which is a measurable 
quantity describing particle distribution at long times after a certain prepared distribution. 
Such diffusion coefficients were measured in Ref. ((J using a fluorescence apparatus. 
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FIG. 4: Graphical representation of recursion relations for the diffusion propagator Q which de- 
termines diffusion coefficient D and vertex proportional to A. Single plain lines are propagators 
G, open circles are averaged force correlators (jSJ) on the momentum shell Ae~ l < q < A, and bare 
vertices are given by —iXki according to the last term in Eq. (|15jl . Single wiggly lines are diffusion 
propagators Q with momentum cutoff A and double wiggly lines are diffusion propagators with 
scaled cutoff, to order A 2 . Dashed circle is the renormalized vertex, to order A 3 . The intermediate 
momenta are integrated over the shell Ae~ l < q < A. 

We must now solve Eq. (J3"4"J) using v r (k) given by Eq. (J3Tj) as an input (and let us not 
approximate v r « v for the moment). This is in general a formidable task analytically. 
When a ~ 0, however, equations simplify considerably, and we get for the small-A; limit 
(assuming a A ~ 1) 



A J — In (AS r 



57T 



which is the same as Eq. (114)1 apart from the numeric prefactor A 



(35) 



17 - l)/2, in 

agreement with Ref. llj. Using the k = limit of Eqs. (fl"4"|) and (|33|) . we find that v r and 
D T essentially satisfy the relation up to a numeric factor of 2/ A w 1.3, if we set I* = 5 r 
and assume aA ~ 1. At low temperatures, however, u r ~ v and v r ^> D r [as the diffusion 
coefficient is in general proportional to T according to the Einstein relation (J2J)]. Solving 
Eq. d3H) then gives Eq. (|T7)jl for 0. 
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